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We investigate the long-time evolution of weakly perturbed single-site breathers (localized sta- 
tionary states) in the discrete nonlinear Schrodinger (DNLS) equation. The perturbations we con- 
sider correspond to time-periodic solutions of the linearized equations around the breather, and can 
be either (i) spatially localized, or (ii) spatially extended. For case (i), which corresponds to the 
excitation of an internal mode of the breather, we find that the nonlinear interaction between the 
breather and its internal mode always leads to a slow growth of the breather amplitude and fre- 
quency. In case (ii), corresponding to interaction between the breather and a standing- wave phonon, 
the breather will grow provided that the wave vector of the phonon is such that the generation of 
radiating higher-harmonics at the breather is possible. In other cases, breather decay is observed. 
. This condition yields a limit value for the breather frequency above which no further growth is 

possible. We also discuss another mechanism for breather growth and destruction which becomes 
P_i ' important when the amplitude of the perturbation is non-negligible, and which originates from the 

oscillatory instabilities of the nonlinear standing-wave phonons. 

I. INTRODUCTION 

Cu 

The concept of nonlinear self-localization is of importance for many physical phenomena, and has appeared in 
a number of different contexts since the pioneering work by Landau Q on the polaron problem in the 1930s. In 
recent years, much attention has been devoted to studies of spatially localized and time-periodic vibrational modes in 
anharmonic lattices (see e.g. || for recent reviews). The general existence of such modes, which have been termed 
discrete breathers, or intrinsic localized modes, as robust solutions to nonlinear (and in general non-integrable) lattice- 
equations was suggested in 1988 by Takcno et al pj. Later, their existence was rigorously proven under rather general 
conditions by MacKay and Aubry j^j by considering the limit of uncoupled oscillators (the so called anticontinuous or 
anti-integrable limit). By means of the implicit function theorem, they showed that the trivial solution of a single-site 
localized vibration at the uncoupled limit could be continued into a localized breather solution for non-zero coupling 
^ between the oscillators, provided that the individual oscillators are anharmonic, and that no multiples of the breather 
frequency resonate with the bands of linear excitations (phonons). As was demonstrated first in jfj), the ideas of the 
C^) • rigorous proof can be turned into an efficient numerical scheme to calculate breather solutions to any desired accuracy. 

Since the discrete breathers appear under very general conditions in anharmonic lattices and provide efficient means of 
energy localization, they have been proposed as candidates to explain experimentally observed localization of energy 
in many different physical areas, e.g. DNA dynamics 0. 

Although, from a fundamental and mathematical viewpoint, the existence theorems for discrete breathers provide 
an important cornerstone for understanding the dynamics of anharmonic lattices, it is probably of even larger phys- 
ical importance to understand the behaviour of a system close to an exact breather solution. By linearizing the 
lattice-equations around the exact solution, one can obtain an approximate description of the dynamics of weakly 
perturbed breathers, and in particular the linear stability properties determining whether small perturbations will 
grow exponentially or not. It was shown in ||, Q that the simplest, single-site, breathers are generally linearly 
stable close to the uncoupled limit, and numerical investigations using standard Floquet analysis (see e.g. Q) have 
shown that linearly stable breathers typically exist also for rather large values of the inter-site coupling. However, 
when considering time-scales large compared to the breather period, the mere linear stability of a breather does no 
longer guarantee the eternal existence of the breather in the presence of small perturbations, and there are still many 
questions remaining concerning the different mechanisms by which breathers may grow or decay, or possibly finally 
be destroyed. If the breathers have a finite life-time, the determination of this life-time is of large importance for 
understanding the role of breathers in real systems. 

It is the purpose of this paper to investigate in more detail some mechanisms for breather growth and decay in a 
simple model system, the discrete nonlinear Schrodinger (DNLS) equation. The DNLS equation is generic in the sense 
that it describes slowly (in time) varying modulational waves in discrete systems in a 'rotating-wave' approximation 
(see e.g. @JT3|); however, due to its extra symmetry properties (see Sec. |TJ) it exhibits some nongeneric features 
among discrete systems, e.g. exact quasiperiodic breathers jy]. The single-site breathers of the DNLS equation 
are stationary states which are linearly stable for all inter-site coupling, and which reduce to the NLS soliton in the 
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continuum limit (see e.g. |5|,|Il|-|l5fl ) . An important application appears in nonlinear optics, where the single-site DNLS- 
breather describes a discrete spatial soliton in an array of weakly coupled waveguides []l6l[l7|| ; recent experimental 
observations Q confirm the successful use of the DNLS-model in this context. 

Some recent numerical investigations jl9) have shown that DNLS-breathers can be spontaneously created from 
noisy backgrounds, in a similar manner as was previously observed for Klein-Gordon (n],|2^] and FPU lattices. 
Typically, this spontaneous energy localization was observed to occur in two steps. In the first step, a large number of 
small breathers are created as a result of the modulational instability |^,|l^,^2|,^3| of travelling plane waves occurring 
for certain wave number regimes. The second step proceeds by inelastic collisions between the breathers, in which 
systematically the big breathers grow at expense of the smaller ones. Thus, the outcome will be a small number 
of large breathers, together with some remaining background of small-amplitude (phonon) oscillations. However, it 
can generally not be concluded from the numerical simulations that this is the true final state of the system, and 
actually long-time simulations for FPU-chains pl| revealed also a third step, in which the interaction with the phonon 
oscillations leads to the final destruction of the breather and equipartition of energy. Thus, to elucidate the nature of 
the final states for typical initial conditions in anharmonic chains, it is necessary to obtain a better understanding of 
the mechanisms for interactions between breathers and small-amplitude perturbations. 

In this paper, we take the following approach. As initial state, we consider an exact single-site breather solution, 
and add a small perturbation corresponding to a time-periodic eigensolution to the equations of motion linearized 
around the breather. These solutions, which can be either localized or extended in space, constitute a complete set 
in which an arbitrary initial perturbation can be expanded. The localized solutions correspond to internal modes 
of the breather Jl3| , p4| -[2^] , while the excitation of an extended solution corresponds to a standing- wave (i.e., non- 
propagating) phonon interacting with the breather. In Sec. O we describe the model and outline the perturbational 



approach which forms the analytical backbone for the interpretation of the numerical results presented in Sees. Ill 



and IV. Sec 



while Sec. IV 



[II discusses the long-time consequences of the interaction between the breather and its internal modes, 
concerns the interaction between the breather and small-amplitude standing-wave phonons of different 
wave vectors. We will find that in both cases, a simple argument based on the conservation laws can be used to obtain 
a sufficient condition for breather growth. In Sec. [VJ we discuss another type of mechanism for breather growth and 
destruction which becomes appreciable when the amplitude of the standing wave is non-negligible (and consequently 
the perturbational approach can be expected to fail), and which has its origin in the recently discovered oscillatory 
instabilities of the nonlinear standing- wave phonon themselves [B7j . Finally, we make some concluding remarks in 



Sec. VI 



Concerning the numerical simulations of the dynamics presented in this paper, unless otherwise stated they always 
apply for a system of infinite size (finite size systems are considered only in Sec. 0). The simulations have been 
performed either by using very large system sizes or by appending damping regions of various sizes to the boundaries; 
in all cases we have carefully checked that the boundary conditions have no essential influence on our results. 

II. MODEL AND FRAMEWORK FOR THE PERTURBATIONAL APPROACH 

A. Model 

We consider the following form of the DNLS Hamiltonian with canonical conjugated variables {iip n }, {"0n} : 

^(Wr l },W; i })=E( C l^+l-^| 2 -^IV^| 4 ) ^E^n- (1) 

rl n 

This yields the DNLS equation 

dl~L 

tyn = TT— = -C(lp n+ l + 1p n -i - 2lp n ) - \lp„\ 2 1pn, (2) 

which, in addition to the Hamiltonian (|l|), also conserves the total excitation norm (or power in nonlinear optics 
applications), 

A^^hAnl 2 ^^"- ( 3 ) 

n n 

The conservation laws for the norm and Hamiltonian are, through Noether's theorem, related to the invariance of the 
DNLS equation (or, more precisely, of its corresponding action integral) under infinitesimal transformations in phase 
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(VVi — ► ipne 1 *) an d time (t — ► t + e), respectively. Defining the 'norm density' A/"„ and 'Hamiltonian density' Ti n as in 
Eqs. (0) resp. (|l|), the conservation laws can be expressed in terms of continuity equations as 

^ + (Jat)„ - (J^) n _ x = 0, (4) 



with the (norm) current density 



and the Hamiltonian flux density 



^ + (^) n -(J«)„_ 1 = 0, (5) 



JaA = 2CIm (6) 



Jn = -2CRe [Vn+itC+i - C)J . (7) 

respectively. These conservation laws are discrete analogs to those existing for the continuous NLS equations with 
general nonlinearities (see e.g. p8| ), however there is no discrete counterpart to the momentum conservation law since 
the discrete equation has no continuous translational symmetry in space. Furthermore, we note that the transformation 
C — ► — C in (|2|) is equivalent to ip n — > (—l) n e~ l4Ct 7p n , and thus we will for the rest of this paper only consider C > 
without loss of generality. 

The single-site DNLS-breather is a stationary-state solution to Eq. (||) of the form 

il> n (t) = 0n(A)e lAt , (8) 

where the time-independent shape {4> n } depends on the frequency A and is spatially localized with a single maximum 
at a lattice site. The breather exists for all A/C > 0; the limit C — > (or A — > oo ) corresponding to the anticontinuous 
limit, where {4> n } is localized at a single lattice-site, while the limit A/C — > corresponds to the continuous limit, 
where {<j> n } approaches the NLS soliton. The single-site breather is a ground state solution to Eq. (||) in the sense 
that it minimizes the Hamiltonian ([[]) for a fixed value of the norm (||), i.e., STL + ASAf = 0, where the frequency 
A appears as the Lagrange multiplier (see e.g. fLqj ). The norm for the single-site breather, Af$, is known to be a 
monotonously increasing function of A, while the Hamiltonian, H.^, is negative and monotonously decreasing (see e.g. 
|p9| , j30f ). From the minimization condition, these functions will be related as 

~dA ~ A ^T (9) 
To describe the dynamics close to the breather (^), we introduce the following perturbation expansion: 

Mt) = {<Pn + Mi) + x2 Vn(t) + X 3 Ut) + A 4 0„M + ...} e l I Adt , (10) 

where e„(0) is the initial perturbation and rj n (0) = £ n (0) = 9 n (0) = ... = 0. Thus, as for the usual stability analysis 
of stationary states (see e.g. ]l^ , ^l| ), the perturbation is applied in a frame rotating with the breather frequency A. 
Substituting into Eq. (^) and identifying coefficients for consecutive powers of the small parameter A yields an infinite 
set of equations, which from 0th to 4th order read: 

- A(j) n + C{(j) n+ i + (j) n -i - 2(j> n ) + \4> n \ 2 <j)n = (11) 

£(A) • {e„} = {ii n + C(e n+1 + e n _i - 2e n ) + 2|0„| 2 e„ + <fae* - Ae n } = (12) 

£(A).{ry„} = -Ce ? 2 1 -20„|e„| 2 (13) 

£(A) • {£n} = -2<P* n tnVn - 20„ (e* n n n + e n T]* n ) - |e„| 2 e„ (14) 

£(A) • {0 n } = -24> n (e„e: + e*Jn + \Vn\ 2 ) ~ 4>l {^n + rfc) - e 2 n r,* n - 2|e„| 2 r; n , (15) 

where the operator C(A) (which is linear over the field of real numbers) is defined from the first equality in Eq. (p^). 
The 0th order equation (O) gives the breather shape {<j) n } (which for the single-site breather can be assumed real 
and positive without loss of generality), while the 1st order equation ( |l2| ) is the linearization of the DNLS equation 
around the breather. 
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B. Solutions to the linearized equations 



To obtain the solutions to the linearized equations ( |l2| ) , we proceed in a similar way as is usually done for continuous 
generalized NLS models (see e.g. |f32"|-j34[ ) and introduce a substitution of the form 



e n (t) ee e«(t) + ie<£\t) = \a (U n + W n ) 4 + \a* ([/* - W*) (16) 

so that 

6M(t)=Re(e„(i))=Re(aC/„e-^ t ) , e$(t) = Im (e„(t)) = Im (aW^-^) . (17) 

Substituting ( |l6| ) into Eq. ( |l2] ) and assuming </>„ real yields 

£ Q W n ee -C(W„+i + W n -i - 2W„) - 4> 2 n W n + AW n = uo p U n (18) 

&U n ee -C(£/ n+1 + U n -x - 2U n ) - 34> 2 n U n + AU n = to p W n , (19) 

where the operators Cq and C\ are Hermitian. Thus, we can obtain the eigenfrequencies lj p and the corresponding 
eigenvectors {{U n }, {W n }) from matrix diagonalization, 

M (0) A {Un}\ = (0 C \( {U n } \ / {[/„} \ 

I {Wn} J ~ \ Ci Oil {W n } J ~ M {W n } J ■ [M) 



To make connection to earlier work MM , we remark that the vector 



{en } }, {<$}) = ({U n }, {-iW n }) is an eigenvector 



of the Floquet matrix with eigenvalue e luJpT , where the time-period T here is arbitrary since the operators Co and 
C\ are time-independent. (The symplectic Floquet matrix is e MpT , where Mj? is obtained from M^ ' by changing C\ 
into — Cx.) Thus, ( |l6| ) is the linear combination of two complex conjugated Floquet eigensolutions which makes ei^ 
and £n real. For the single-site breather, all eigenvalues lu p of M( ' are always real, implying the linear stability of 



the breather for all parameter values A/C > |35|. Accordingly, we can also choose the eigenvectors ({U n }, {W n }) of 
M^ ) to be real and normalized, in which case the phase of the amplitude a describes the symmetry properties of the 
solution ( |l6|) under time reversal: choosing a real yields a time-symmetric solution, e n (t) = e* (— t), while choosing a 
purely imaginary yields a time-antisymmetric solution, e n (t) = — e*(— t). 

For an infinite system, the spectrum of the (non-Hermitian) matrix M^ ' can generally be divided into a continuous 
(phonon) part, corresponding to extended eigenvectors, and a point spectrum corresponding to localized eigenvectors. 
The phonon spectrum for any localized solution {</>„} is easily obtained from the limit |n| — > oo, since the condition 
4> n — > implies that the operators Cq and C\ become identical and Eqs. ( |l8| ) - ( |l9| ) reduce into two uncoupled 
equations for the linear combinations a„ = U n + W n resp. b n = U„ — W„- Assuming a n ~ e ±ll?a ™ and b n ~ e ±lqbn , 
respectively, yields the dispersion relations 

uj p = A - 2C(cosq a - 1), (21) 
uj p = -A + 2C(cosg b - 1), (22) 

from Eqs. ( [l8| ) - (|l9|). Thus, the continuous spectrum of the matrix M^ ) consists of two branches, symmetrically 
located around ui p — 0, and since A > for the single-site breather these two branches never overlap. Note also 
that two eigenvectors with eigenvalues ±oj p correspond to the same solution to Eq. ( |l2] ) (changing the sign of ui p in 
Eq. ( |l6|) is equivalent to changing U n <-> U*, W n ^ —W*, a <-> a*), and therefore it is enough to consider e.g. ui p > 0, 
in which case b n — U n — W n always vanishes exponentially as n — > ±oo. 

When A/C is not too large, the linear spectrum around the single-site breather contains also two pairs of nonzero 
isolated eigenvalues uj p , which correspond to the two internal modes of the breather [ ^3p^ |. One of these modes is a 
spatially symmetric, 'breathing', mode, while the other is a spatially antisymmetric, 'translationaP or 'pinning' mode. 
Numerically, it has been found that the breathing mode exists for < A/C < 1.7, while the pinning mode exists for 
< A/C < 1.1. The numerically calculated internal mode frequencies uj p as a function of breather frequency A are 
shown in Fig. [l]. Note that as A/C — > 0, the breathing mode frequency approaches the lower edge of the phonon 
band (but always stays outside the band (26)), while the pinning mode frequency approaches zero (but always stays 
nonzero). This is consistent with the fact that the soliton solution of the continuous NLS equation has no breathing 
mode (due to its exact integrability) , and has a translational mode with zero frequency due to the translational 
symmetry of the NLS equation. 
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To obtain a complete set of solutions to Eq. ( pL 2[ ) in which an arbitrary initial perturbation e n (0) can be expanded, we 
must include also the zero-frequency solutions, which generally can be written as a superposition of two fundamental 



modes. One of these modes ('phase mode' j36| ) is the solution W n = <fr n to the homogeneous equation (|18|), CoW n = 0. 
The corresponding perturbation e„ = i(j) n describes a rotation of the overall phase of the breather. The second mode 
('growth mode' |36|) is obtained by solving the inhomogeneous equation C\U n — —4> n , which has the solution 
U n = d(j} n /dA. The corresponding solution to Eq. ([l2] ) is e„ = d(j) n /dA + i(j) n t, and corresponds to a time-linear 
growth of the perturbation representinga small change in the breather frequency. 

Although the set of eigensolutions ( |l6| ) together with the two zero-frequency modes forms a basis for the space 
of solutions to Eq. ([l2]) (there are no bifurcations, which could result in additional 'marginal modes' |3(| with time- 
linear growth at degenerate eigenvalues), this basis is in general not orthogonal using the ordinary scalar product, 
and typically there is a considerable overlap between the solution corresponding to the internal breathing mode and 
the zero-frequency modes. However, in analogy with e.g. J32|-Q, we can define a 'pseudo-scalar' product between 

any two vectors {V^ 1} }) and ({C/^}, {wl 2) }) by 

^(f/(%( 2 )*+lfW[/f). (23) 

n 

This product is formally not a true scalar product, since for the general case the product of a vector with itself 
as defined by (P3|) is not necessarily positive. However, when ({U n },{W n }) is a real eigenvector of M^ ) we have 
E>> ^n^n = tj-^2 n W n £oW n from ([l^), and the operator £ is positive definite for all W n ^ <p n Q. With this 
product, it follows from ([l8]) - ( p^ ) that all eigensolutions with different (real) eigenfrequencies lo p are 'orthogonal' in 
the sense that 

( W W - £ (t/Wwf >* + W^r) = 0, (24) 

n 

and the only nonzero product involving the zero-frequency modes is the cross-product between the phase mode and 
the growth mode (3f| : 

— ^" ld ^>0, (25) 



<9A 2 dA 



where is the norm (||) of the breather with frequency A. We also remark that, since the product ( p3|) multiplied 
by a factor i is just the symplectic product between the two vectors ^{4^^}, {cn^^H and ^{el r ^ 2 ^}, {en^H , 
the sign of the product of an eigenvector with itself can be interpreted as the negative of the Krein signature of the 
corresponding pair of Floquet eigenvalues H . 



C. Strategy for the perturbational approach 

We now consider as initial state a single-site breather perturbed in the direction of a single eigenmode (|l^) (lo- 
calized or extended) of the linearized equations fll^), and wish to describe qualitatively the long-time effects of this 
perturbation using the expansion (|l^). In general, taking into account terms up to order p in this expansion yields 
a solution to the DNLS equation which is correct to 0(X P+1 ), i.e., for long but finite time-scales for small initial 
perturbations (note that the expansion parameter A plays the same role as the mode amplitude a.). As is wellknown 
however, this kind of expansion in general diverges due to resonances between solutions to the homogeneous equa- 
tion (|l2l) and the inhomogeneous terms appearing in the right-hand sides of Eqs. (|l3|)-([l5|) and the corresponding 
higher-order equations. A resonance with a solution belonging to the continuous spectrum results in a bounded but 
nonlocalizcd solution corresponding to outgoing radiation, while a resonance with an eigenfunction belonging to the 
discrete spectrum gives a spatially localized response which diverges linearly in time. However, up to any finite order 
these divergences can be systematically removed by allowing a slow time dependence of the independent variables, 
which in our case are taken to be the mode amplitude a and the breather frequency A. This procedure adds additional 
terms to the equations, which can be tuned so that the divergent parts of the response disappear. In other words, 
these two quantities are used as collective variables which, together with the outgoing radiation fields, are expected 
to describe the main features of the asymptotic dynamics if the initial perturbation is sufficiently small. 

The second order correction is given by the inhomogeneous Eq. (|l3|), which with the substitution ( f[6| ) becomes 
(choosing tj> n , U n , and W n real without loss of generality): 
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£(A) • { Vn } = -~4> n [\a\ 2 {W 2 n + Wl) + {Wl Wl) Re (aV 2 ^*) 



-2iU n W n Im (a 2 e 



-2lUJr>t 



)]• 



(26) 



Thus, the right-hand side contains one static part and one part involving the frequencies ±2w p . It acts as a periodic 
force with frequencies and 2tu p , and since all terms contain the factor </>„ this force is localized at the breather region. 
The response to this force will remain bounded and localized unless the corresponding homogeneous equation ( |l2| ) has 
a solution with frequency Oresp. 2ui p which is non-orthogonal to the corresponding part of the right-hand side in d2 



As will be shown in Sec. Ill B| , a non-zero overlap between the static part of ( p6| ) and the zero-frequency solutions of 
( |l2| ) yields a (time-independent) shift of the breather frequency. Moreover, if A < 2\u p \ < A + 4C, so that 2uj p is inside 
the phonon band of the homogeneous equation, a resonance will generally occur, resulting in radiation with frequency 
2uj p emitted from the breather region. The strength of the radiation field is determined by the (generally nonzero) 
overlap between the 2a; p -part of ( |26| ) and the corresponding homogeneous solution (see Sec. IIIB). In a similar way, 
we obtain that the right-hand-side of the third order equation ( |l4| ) contains the frequencies u> p and 3w p , the fourth 
order equation ( [l5| ) contains the frequencies 0, 2ui pi and 4w p , and in general the pth order equation contains as its 
highest harmonic the frequency puj p . Accordingly, we conclude that if 



A < p|w p | < A + 4C, 



(27) 



so that ptOp belongs to the phonon band, the perturbed breather will radiate to pth order. The consequences of 
this radiation for the breather itself will be discussed in Sees. Ill and IV for the cases of localized and extended 
perturbations {e„}, respectively. 



III. BREATHER INTERACTING WITH INTERNAL MODES 



With the initial perturbation e„(0) of the single-site breather corr esponding to a spatially localized eigenmode of the 
linearized equations (^|), it is clear from the discussion in Sec. II C that higher order radiation always will be created, 
since for any internal mode frequency lo p there is always an integer p such that pw p belongs to the phonon band 
and ( ^7|) is fulfilled. Moreover, from the numerical results presented in Fig. [l| we find that the spatially symmetric 
breathing mode always radiates to second order, since (|27]) always is fulfilled for p — 2, while the antisymmetric 
pinning mode radiates to second order only when A/C > 0.480. Thus, due to this radiation from the breather the 
total norm contained in any finite region around the breather(i.e., the total norm of breather -I- internal mode) will 
always decrease with time. However, the main concern here is the long-time effect of the internal-mode excitation 
on the breather itself, and thus we must investigat e whe ther there will also be some transfer of energy between the 
breather and its internal mode. We will first (Sec. [II A) show results from direct numerical integration of Eq. (0); 
then we willgive two alte rnativ e approaches to the analytical interpretati on of these results based on the higher-order 
equations ( 13 )-(p~5[) (Sec. Ill B| ) resp. the conservation laws (||)-(||) (Sec. IIIC). 



A. Numerical simulations 



In Fig. we show a typical example on the long-time evolution of a breather when the initial perturbation is taken 
in the direction of its internal breathing mode. As is seen from Fig. ^ (a), the amplitude of the breathing mode decays 
slowly with time as a consequence of the losses due to generation of second order radiation, and a careful study of its 
envelope \a(t)\ indicates that it decays as 



\a(t)\ 



Ko)| 



Vi + 7Ko)| 2 * 



(28) 



where 7 > is a constant. This is consistent with a similar result obtained for the continuum NLS equation with 
generalized (non-cubic) nonlinearity 33 1; the analytical motivation for this result (which is analogous to that of the 
continuum model given in Ref. p3|]) is given in the following subsections. 

However, the main result of this section is illustrated in Fig. || (b) and (c). Fig. ^ (b) is obtained by calculating the 
time-average of the central-site intensity as 



(l^o! 2 ) t 



^I>n (£ fe )| 2 , 
fe=l 



(29) 
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where tu is a set of closely spaced time-instants. It is clear that the interaction between the breather and its internal 
mode asymptotically leads to an increase of the average peak intensity, i.e., to breather growth. The same phenomenon 
is illustrated also in Fig. || (c), where we have plotted the difference between the instantaneous breather frequency 
calculated at time t, A(t), and the frequency of the unperturbed breather Aq. From this figure, we can also conclude 
that there are two different mechanisms causing the shift of breather frequency. Firstly, there is an initial (almost 
instantaneous) rather large frequency shift, which can be interpreted as an adaption of th e initi ally perturbed breather 



to the breather which is 'closest' to the initial condition. As is shown below in Sec. 1IIB, this time- independent 



frequency shift, which is observed to be always positive for the breathing mode, is a consequence of the overlap 
between the static part of the right-hand side of the second-order equation (^6|) and the zero-frequency modes. 
Secondly, there is the slow, continuous increase of the breather frequency which corresponds to the slow increase of 
(\ipno | 2 ) t m Fig- § (b), indicating a continuous transfer of norm from the internal mode to the breather. It is described 
by the static part of the right-hand side of the fourth-order equation (|l5|) (see Sec. Ill B| ). 

When the initial perturbation of the breather is taken in the direction of its internal pinning mode we observe, 
just as for the breathing mode, that the breather-internal mode interaction asymptotically always leads to breather 
growth. An example is shown in Fig. where the parameter values have been chosen so that the lowest harmonic 
that enters the phonon band is 3w p (A/C = 0.45 < 0.480), We observe two qualitative differences compared to 
the case with breathing mode excitation. Firstly, since in this case the first phonon resonance occurs only in the 
third order equation (J14|) , the decay of the internal mode amplitude will be slower, and a good fit is obtained by 

\a(t)\ ~ |o(0)| (l + 7 |a(0)| 4 i)~ 1/4 . This agrees with the general result when puj v is the lowest harmonic that enters 
the phonon band obtained for the continuum generalized NLS e quation in Ref . |33| ; the derivation of the corresponding 
result for the discrete case (see Eq. (ffTj)) is given in Sec. Ill C . As a consequence of the slower decay of the internal 
mode amplitude, the breather growth will also be slower when p > 2, as can be seen from Fig. |^ (c) by comparing the 
time-scales with those of Fig. g. Secondly, the initial shift of the breather frequency for a pinning mode excitation 
is always much smaller than for the breathing mode excitation (also when 2lo p is i n the phonon band), and is when 
A/C > 0.55 also observed to be negative. The explanation for this is given in Sec. Ill B . However, it is important to 
stress that also in the cases where the initial frequency shift is negative, we find that the continuous breather growth 
always will give an asymptotic frequency shift which is positive. 



B. Analysis of higher-order equations 

Here, we will analyze the higher-order equations (P"3"l)- ([l5|) by making use of the strategy of systematically removing 
the appearing divergent parts as outlined in Sec. fl C\ (in analogy with the treatment of the continuous NLS-type 
equations in e.g. Refs. J32|-p4]|). First, we show how the dominating contribution to the time-independent frequency 
shift observed in the numerical simulations above can be calculated from the static part of the right-hand side of the 



second-order equation (13). This frequency shift can be explicitly taken into account by replacing A in Eq. (|10|) with 
Aq + A 2 A2, where Ao is the unperturbed breather frequency and A2 the second-order shift to be determined. This 
implies that the additional term A2(j> n will be added to the right-hand side of Eq. (p6|). Writing the response to the 
static part of (p6|) as ?/i s ' = ja! 2 ^! 8 - 1 + iwn ^) with real Un and then yields 





with M' ' as defined by Eq. (pcf). If the expansion of the right-hand side of ( |30| ) in the complete set of vectors 
consisting of the eigenvectors of M^ ) (including the phase mode) and the growth mode contains some component 

(s) 

on either of the two zero-frequency modes, the response rfa will not remain bounded but diverge linearly with time. 
Thus, in order to remove this divergency, the frequency shift A2 must be chosen so that both these components are 
identically zero. The component on the growth mode is trivially zero, while the component on the phase mode is 
obtained by applying the pseudo-scalar product ( p^ ) with the vector corresponding to the growth mode and using 
(p5|). Demanding this component to be zero yields 



dA n 



This is typically positive for the breathing mode, since the dominating contribution to the sum in ( pf|) comes from the 
central site no, and is always positive. For the spatially antisymmetric pinning mode, there is no contribution to 
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this sum from the central site, since U na and W no are zero. The change from a positive to a negative frequency shift 
when increasing A in this case is related to a qualitative change of the nature of the growth mode 5^-: for A > 0.55 
we find that < for all n ^ no, so that all terms in the sum in Eq. d3l]) are negative, while becomes positive 



g\ \ yj iui ein it -f- ou unc^L cm ijcimo ill Luc ouiii ill V|^J/ negative, vviiiic 

also for sites in the neighborhood of uq for smaller A. 

For the rest of the analysis in this subsection, we assume for calculational simplicity that the internal mode fre- 
quency is such that 2uj p is inside the phonon band (and thus it is not applicable for the pinning mode excitation 
when A/C < 0.480). Then, the non-static part of the right-hand side of the second order equation (|26| ) will gen- 

-2lLJ v t i 



erally give rise to a non-localized response, which can be written in the form = \a 2 + ^ e 

■ia* 2 (un^ — WrP e 2j "p*. This response corresponds to the radiation field going out from the breather region, and 

since the right-hand side of ( |26| ) is spatially localized and symmetric, this field should asymptotically correspond to 
two identical linear waves propagating to the left (right) for n — > — oo (+oo). Thus, the boundary conditions can be 
written as 



(2) , w ™ -> r 2 e ±lq2n , n ±oo, (32) 



with q2 = arccos ^1 — —§cj according to (pl|). Defining for general u> the matrix (cf. (pp|)) as 

the functions zti 2 '* and wi 2 -* are seen from ( |2(f ) to be determined by 



MH , _^ , (33) 



(34) 



M (2. p ) / {ttf'} ) _ <f>n ( {2U n W n } \ 



Since for general u> every eigenvector of with eigenvalue /i is also an eigenvector of M^ ^ with eigenvalue (x + uj, 

the right-hand side can be expanded on the basis of eigenvectors of M^ ^ (including the zero- frequency modes). The 
strength of the radiation field is then given by the expansion coefficient for the (continuous spectrum) eigenvector of 
M^ ) with eigenvalue 2w p , since this corresponds to the eigenvalue /i = of M' 2 " p ', and thus a spatially non-bounded 



response in (34). Using the orthogonality relation (|24j), this coefficient is simply the 'overlap' between the right-hand 
side of ([34|) and the eigenvector of M( ) with eigenvalue 2uj p calculated with the pseudo-scalar product (|2^), which is 
generally nonzero. 

Next, we show how the dominating contribution to the decay of the internal mode amplitude as given by Eq. ( pq ) 
is obtained from the condition that £ n in the third-order equation ( |l4| ) should remain bounded. To this end, we 
assume a slow time-dependence of the internal mode amplitude of the form a = a(A 2 t), and consider the response to 
the terms with frequency uj p in the right-hand side of Eq. ( pT[ ) 'corrected' by the additional terms (— ia + A 2 a)(£/„ + 
W n )e~ ZUJpt + (—ia* + A2a*)([/ n — W n )e IWpt appearing as a consequence of including the time-dependence of a and the 
second-order frequency shift A2 from Eq. ( |3l"| ) in the perturbation expansion @. Writing the response to this part 

as $J = ±\ a \ 2 a (u ( n 3) + w { n 3) ^ e"^ + \\a\ 2 a* (u {3) * - u4 3) *) e^ 4 yields 

,(3h 



i l W~a W n + <t>n (3C/„^ 2) + W n w£ ] + 6C/„^ s) ) + \ (3U, 




(35) 



A bounded response for exists only if the vector on the right-hand side of (|3^) has no component in the direction 
of the internal mode eigenvector ({£/„}, {Wn}), since this is the eigenvector corresponding to the eigenvalue zero of 



M( w p). Using the orthogonality relation (24), this component is obtained by application of the product ( p3|) with 
the vector {{U„}, {W n }), and the condition that this component must be zero determines the time-evolution of a. 
Considering only the absolute value |a| 2 , the resulting equation has the form 4z (\a\ 2 ) + 7H 4 = 0, which has the 



> 0, (36) 



desired solution (|28j). The constant 7 is given by 

E n K [2U n W n lm (wl?) + (3U 2 - W 2 ) Im (u^)] 8CWp | r2 |2 sinq2 



7 = 



E„ u n w n e„ w n c w„ 
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where the second equality is o btained using Eqs. ([T^), (^) and (0), and the positivity of 7 follows from the fact 
that, as mentioned in Sec. II B , the operator £0 is positive definite for all W n ^ <f> n |3^ |. 

Finally, we show how the continuous increase of the breather frequency appears from the divergent response to the 
static part of the right-hand side of the fourth-order equation (fL5h. In its unmodified form, this part is given by 



Ma| 4 [6U n Re (u^ + 2W n Re («;<?>) 
+6(«W) a 



I- 2iI7 n Im 
3|<4, 2) l 2 + k 



,(3) 



(2)|2 



2ilm 



2iW n Im (u^) 

} 



{2 (3(7 2 + W 2 ) «W + (3f/ 2 - W*) Re («£0) 



2(7„VK„Re 



(-i 2) ) 



-2iU n W n lm it 



,(2) 



+ z(C/ 2 -3W 2 )lm 



(37) 



Now, it is clear from ( |3l| ) and ( |28D that the time-dependence of a will induce a time-dependence of the second- 
order frequency shift A2 of the form A 2 (A 2 i), so that we can express the total breather frequency up to order A 4 
as A(t) = Ao + A 2 A 2 (A 2 t) -I- A 4 A4, where a 4th order correction A4 has also been included. Then, we must also 
take into account the time-dependence of the breather shape <p n by writing </>„ (A (<;)). As a consequence, the term 



<9A 



A 2 + A 4 0„ will be added to the expression (pT\) for R 



(is) 



to this total force as 



l«| 4 (1 



(is) 



,.(4«) 



with real and wif s 



in the right-hand side of (|15|), and writing the response 
yields 




— Im 

- Re 



K s) )} 
K s) )} 



(38) 



The response On^ will be bounded in time only if the right-hand side of ( |38| ) has no component either on the growth 
mode or on the phase mode, which gives two conditions for the determination of A 2 and A 4 . Using (p3|)-(p5|), we 
obtain by demanding the expansion coefficient for the growth mode to be zero: 



A 2 = 



J dA n 



W n lm I u 



,(3) 



2U n W n lm 



U n Im I w 



,(2) 



,(3) 



Im 



{ZW 2 n -Ul) Im (™i 2) )]} 



(39) 



(Similarly, A4 is obtained by demanding the component on the phase mode to be zero.) Thus, the dominating 



contribution to the frequency growth should be of order A 2 
( p8| ) of the internal-mode amplitude we obtain qualitatively 

A(f)-A ~|a(0)| 2 (d-C 2 



so that with the approximate time-dependence 



(40) 



l + 7 |a(0)| 2 i 



which is in good agreement with the numerically observed time-dependence of the frequency shift as shown in Fig. |2j 
(c). However, the positivity of A 2 is not easily seen from the expression (^9|), and therefore we will in the following 
subsection derive an alternative expression from which the positivity follows immediately, using the conservation laws 
for the norm resp. Hamiltonian. 



C. Approach using conservation laws 



We consider first the conservation law (^) for the total norm (|^) contained in any large but finite region around 
the breather. Averaging over a time-interval [t,t + 2tt/lo p ] and using Eqs. (|l(i|), (|l6|), we can write the time-averaged 
norm to second order in \a\ (putting A = 1) as 



Wt(*) 



E 



IK + wi 



(41) 



(Note that there will be no contribution at order 2 from the static second-order correction Un \ since the renormal- 



ization of the breather frequency A according to (|3l|) yields J2 r 



0.) In the case when 2w p belongs to the 



phonon band, we obtain the following balance equation, which is correct up to order 
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where we have used Eqs. (|4|), (|6j) and (|32|), and 52 is as defined below Eq. (p2[). 

Similarly, we can use the conservation law (|^) for the Hamiltonian ([l]) together with the general expression for the 
Hamiltonian flux density (]?]) for a small-amplitude plane wave ip n = Ae l ^ qn ~ wl ^ li '\ 

J H = 2\A\ 2 Cu(q)smq, (43) 

to write the balance equation for the total time-averaged Hamiltonian in the same region for the case of second order 
radiation: 

d(H) t dHd, ■ d(H) t d\a\ 2 . . 4 2 

= -^f A + ^ ; ^ = Jw(-oo) - Jw(+oo) = ~4C\a\ i \r 2 \ 2 (2u Jp - A) sing 2 , (44) 

which is also correct to order \a\ 4 . The lowest order contribution to the derivative d g^ji can be obtained using the 
first equality in the equation of motion (||) and its complex conjugate as follows: 

d(H) t _ 1 d(H) t _ / &H_ dip n , dH dtp 



9|a| 2 a* 9a a* ^— ^ \dip n da dip* da 

= -\ A T, K + w ") +^ P J2 u * w " + °(H 2 )' ( 45 ) 

n n 

which is always negative for an internal mode excitation since \u> p \ < A. Thus, we can combine the two balance 
equations ( (42|) and (p4|), and using Eqs. (g) and ( |45| ) we obtain the expression ( |2§| ) for |a(t)| with 7 as in (|36|), 
together with the following expression for the frequency growth rate from which its positivity is immediately seen: 



4CM 4 |r 2 | 2 sin 92 f En K + ™l 

En UnW n 



A = ^ 0111 " 2 ( ^V-n — - 1 ] > 0. (46) 



dA 

This approach also has the advantage that it is easily generalized to the case where the lowest harmonic that enters 
the phonon band is puj p with p > 2, i.e., for the pinning mode excitation when A/C < 0.480. Then, we can write the 
boundary conditions at the infinities to lowest order in a as ip n — ► a p r p e l ^ iiqpn ~^ pulp ~ t ^ t \ n — > ±00, where r p ~ 1 and 
<7 p = arccos ^ A 2 c" P + according to (pl|). Consequently, we can proceed exactly as above, writing down the balance 

equations for the norm and Hamiltonian to order \a\ 2p just by modifying the right-hand sides of Eqs. ([42] ) resp. ( |44| ) 
by replacing |a| 4 with |a| 2p , r 2 with r p , g 2 with q p , and 2w p with pw p . Combining the balance equations yields the 
following general expression for the time-dependence of the internal mode amplitude, 

/.x, _ H0)l _ 4pC|r p | 2 sin q p 

lWl [l + (p-l) 7p \a{0WP-H] 1 ^-V ,1P ~ E n UnW n >U ' (4?j 

which is the analog to the expression obtained with similar arguments in |}3| for the continuum NLS models. And 
most importantly, we obtain a general expression for the breather frequency growth rate which is positive for all p: 

■ AC\a\ 2p \r p \ 2 sing p ( gEj^tg) A n UR) 

A = { 2 ^u n w n - 1 ) >0 - (48) 



Thus, integrating Eq. (|48| ) using the time-dependence (|47| ) of the internal-mode amplitude, we obtain that the domi- 
nating contribution to the frequency growth generally can be written qualitatively as 

A(i) - Ao ~ |a(0)|2 ( Cl - c ^ + {P -ihM0)\^t ) • (49) 

Let us finally point out that the approach used in this section provides a simple physical interpretation of why the 
interaction between the breather and its internal mode in addition to generate radiation also should lead to breather 
growth. Since the expression (f43|) for the Hamiltonian flux density of a small-amplitude plane wave is positive when 
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uj(q) and q have the same sign, the Hamiltonian energy for a plane wave always propagates in the same direction 
as the wave itself. Thus, the second (or higher) order radiation emitted from the breather region will always carry 
away a positive amount of the Hamiltonian energy, or, equivalently, negative Hamiltonian energy will flow into the 
breather region. Moreover, from Eq. jiq ) it is clear that the contribution to the Hamiltonian from the internal mode 
always is negative and a monotonously decreasing function of its amplitude, and thus the decay of the internal mode 
would cause an increase of the Hamiltonian in the breather region. Consequently, since the Hamiltonian of the pure 
breather is a monotonously decreasing function of the breather frequency, the breather should grow so that the total 
Hamiltonian of (breather+internal mode) decreases. A similar mechanism was recently found to cause soliton growth 
in the parametrically driven continuum NLS equation in the regime of oscillatory instability and this type of 
argument has also been used to explain the 'quasi-collapse' of a broad excitation to a narrow localized state in the 
two-dimensional DNLS equation |37j] . 



IV. BREATHER INTERACTING WITH STANDING- WAVE PHONONS 



We choose, as in the previous section, the initial perturbation e„ (0) to be an eigensolution of the linearized equations 
(|l2|), but now we consider the case of a spatially extended perturbation. Choosing without loss of generality a solution 
of the form ( |l6| ) with positive frequency lo p yields the asymptotic behaviour 

U n , W n — * cos(qn ± 8); n — > ±00, (50) 

where the wave vector q (0 < q < 7r) is determined by the dispersion relation (plj), and 8 is the phase shift across 
the breather. Thus, the excitation of an extended eigenmode corresponds to an interaction between the breather 
and a non-propagating, standing-wave phonon with small amplitude a. As was mentioned in the introduction, these 
standing-wave phonons are generally unstable, but since the instabilities become exponentially weak in the small- 
amplitude limit they are expected to have very little effect on the breather for the perturbation sizes and time-scales 
considered in this section. We will return to the effects on the breather of these instabilities in Sec. [v], where larger 
perturbations are considered. 



In contrast to the case of excitation of a localized internal mode discussed in Sec. Ill, where higher order radiation 
always was emitted from the breather region, the condition ( |27| ) yields that for the standing-wave perturbation the 
breather will radiate to higher order only ifA<w p <A/2 + 2C, so that both lo p and 2uj p are inside the phonon band. 
In terms of the phonon wave vector q, this means that there is a critical value q c , 

q c = arccos (Jg) > ( 51 ) 

such that for < q < q c second order radiation will be emitted from the breather region, while for q c < q < tt 
(A/2 + 2C < ujp < A + AC) all multiples of w p are outside the phonon band, and no higher order radiation is emitted. 
As will be shown below, these two regions yield qualitatively different scenarios for the long-time evolution of the 
perturbed breather. We also note that for A > 4C we have q c = 0, so that for the highly localized, high-frequency 
breathers no phonons can generate higher-order radiation (note also that there are no internal modes in this regime). 
Furthermore, we always have q c < 7r/2, so that the regime of higher-order radiation generation is a subset of the 
regime < q < ir/2 where modulational instability for travelling plane waves occurs H. 

Let us first discuss the case q < q c . A typical example of the long-time evolution of a breather interacting with a 
small-amplitude standing- wave phonon with q < q c is illustrated by Fig. As is seen from Fig. ^ (a), the amplitude 
of the oscillations remains essentially constant in time, but a closer inspection reveals that the average value of \ipn 1 2 
asymptotically increases with an apparently constant rate (see inset in Fig. ^ (b)). Similarly, Fig. [| (b) shows that 
also the total norm contained in any finite region around the breather asymptotically increases linearly with time. 
We find that these results are generic for all cases when the phonon wave vector q < q c (the spatial symmetry of the 
phonon is not important for the asymptotic behaviour), and thus we conclude that in this regime, the breather can 
'pump' energy from the phonon (which is infinite for an infinite system), and thereby grow. 

In the same spirit as for the internal mode excitation in Sec. Ill C| , we can give a simple argument based on the 



conservation laws to motivate why the generation of second order radiation should lead to breather growth. To this 
end, we assume that the initial standing-wave phonon is infinitely extended, and that far away from the breather a 
stationary regime will be reached corresponding to the following boundary conditions 
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Thus, we have taken into account the second-order radiation with frequency 2uj p generated at the breather region 
but neglected possible higher-order radiation; moreover the resonance at the original phonon frequency ui p in the 
third-order equation ( |l4| ) has been taken into account by allowing the incoming and outgoing complex amplitudes a 
and r to be different. We can then, in analogy with Eqs. ([l2|) and (Q), write the balance equations for the total norm 
and Hamiltonian contained in a region around the breather averaged over a time-interval [t, t + 2it /oj p \ as 

^dT = ^ A = - ( J *(+°°))t = 4C [(M 2 - M 2 ) sin ? - N* sing 2 ] , (53) 

and 

^ = ^A=(J w (-oo)} t -(JK( + oo)} 4 

= AC [(\a\ 2 - \r\ 2 \)(u; p - A) sing - \r 2 \ 2 (2cu p - A) sin g2 ] , (54) 

respectively. Here, we have used the facts that the time-average of the norm current density (^) and the Hamiltonian 
flux density (Q) are additive quantities for small-amplitude plane waves, and that in the stationary regime, the mode 
amplitudes a, r and r 2 are time-independent. Combining Eqs. ((5^) and (jHJ) and using (^|), we obtain that the breather 
frequency grows with a constant rate given by 

■ 4C|r 2 | 2 sin g2 

A = — — > °- ( 55 ) 

dA 

The physical interpretation of this result is, similarly as for the case of internal mode excitation, that the generation of 
higher-order radiation results in a net flow of negative Hamiltonian energy into the breather region, which is absorbed 
by the breather by increasing its frequency and maximum amplitude. This process is similar to the one observed for 
the two-channel phonon scattering on breathers in a Klein-Gordon model with a Morse potential in Ref. [Q; however 
in the latter case the second outgoing wave resulted from a resonance in the linearized equations and were therefore of 
the same order of magnitude as the incoming wave, and moreover the outcome in this case was breather decay since 
the energy of a Klein-Gordon breather is an increasing function of its amplitude. 

Considering now the regime q > q c where all multiples of lo p are outside the phonon band, the most important 
conclusion from our extensive numerical investigations is that we never observe breather growth. Instead, we sometimes 
(but not always) observe a very slow decrease of (\ipn \ 2 ) t , and an increase of the fluctuations around this mean- value. 
This behaviour is illustrated by Fig. [| (a) and (b). A possible interpretation of these results is that, since the higher 
harmonics which are created by the breather-phonon interaction cannot propagate, they stay trapped around the 
breather. Thus, this could lead to a transfer of energy from the 'pure' breather, which acquires more and more 
internal frequencies and becomes a 'chaotic breather' pj] |. Another possible interpretation is that the increase of the 
oscillation amplitude is connected with the oscillatory instabilities of the standing waves; as we will show in the next 
section these instabilities provide a mechanism for breather decay. However in some cases, illustrated by Fig. |^ (c) 
and (d), the oscillation amplitude as well as its average value apparently approaches a constant limit value. We have 
at present no explanation for this behaviour (as will be discussed in the next section, there exist exact 'phonobreather' 
solutions which could be candidates for such a final state, but they are unstable); it is possible that the time ranges 
that we were able to study with sufficient numerical accuracy in these cases simply were too short to observe the 
scenario described by Fig. || (a) and (b) . 

To conclude this section, we repeat our main result that breather growth is observed if and only if q < q c , where 
q c is given by (|5l]). The fact that q c = for A > 4C thus implies the existence of an upper limit beyond which 
the breather cannot grow with the type of perturbations considered here. We would also like to relate our results to 
recent numerical simulations of breathers interacting with propagating phonons in Klein-Gordon j3{| and FPU pl| ] 
lattices. For the Klein-Gordon lattice with a (soft) Morse on-site potential, phonons with small wave vector q were 
observed to yield breather growth, while phonons with large q caused breather decay. For a FPU lattice with hard 
anharmonicity, the opposite situation was observed, i.e., small-q phonons caused breather decay and large-g phonons 
breather growth. The fact that the situation for the hard FPU lattice was opposite to that of the soft Klein-Gordon 
lattice could be expected, since in the former case the modulational instability occurs for large q, whereas soft Klein- 
Gordon and DNLS lattices with C > are modulationally unstable for small q. However, we stress that the relation 
between plane-wave modulational instability and breather growth is nontrivial, and at least for the case considered 
in this paper the critical value q c for breather growth from interaction with standing-wave phonons differs from the 
critical value q = it/2 for modulational instability of travelling waves. 
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V. BREATHER GROWTH AND DESTRUCTION FROM STANDING- WAVE INSTABILITIES 



In this section, with the aim at describing the interaction between a breather and a standing-wave phonon with 
non-negligible amplitude, we will take a slightly different point of view than in the preceding sections. Instead of 
choosing as initial condition an exact breather solution and adding a perturbation corresponding to an eigenmode of 
the linearized equations, we will here consider initial conditions which are exact phonobreather Q (or nanopteron) 
solutions. By definition, a phonobreather consists of a spatially localized breather on top of a spatially extended tail 
which is a nonlinear, standing-wave phonon. (There are also solutions where the tail is a propagating wave |3J 
we will not discuss them further here.) Phonobreathers exist generically for nonlinear lattice-equations (see e.g. 
but their existence normally requires an integer relationship between the breather and phonon frequencies. However, 
for the DNLS equation phonobreathers exist for any (rational or irrational) relation between the two frequencies [[ll] , 
as a consequence of the additional invariance of the equation under global phase transformations. 

Since the phonobreathers are exact solutions consisting of a breather part and a standing- wave part, one could 




expect them to be attractors for the initial conditions considered in Sec. IV. However, as was shown recently j27| , 
generically for soft Klein-Gordon and DNLS models with C > all phonobreathers with phonon wave vector q ^ it 
will be linearly unstable if, for fixed phonon amplitude a, the linear coupling C is larger than some threshold value 
C cr {a,q) (i.e., away from the anticontinuous limit). (For lattices with hard potentials and DNLS with C < 0, the 
stable phonobreather has q = 0.) These instabilities are caused by an oscillatory instability of the standing- wave 
phonon itself, which can be understood by considering the construction of a nonlinear standing wave with wave vector 
q close to tt at the anticontinuous limit C — by introducing a periodic array of discommensurations or 'defects' in 
the nonlinear phonon with wave vector tt and amplitude a, ip n = o(— l)"e~ 4 ' 4C ~' a )*, which is linearly stable for all a 
and C > ||. In the anticontinuous limit C = 0, each defect consists when n/2 < q < tt of one extra site with tp n = 
added to the 7r-phonon, which consequently suffers an additional phase shift of tt across each defect. For < q < n/2, 
each defect consists of several consecutive zero-amplitude sites with associated phase shifts (a general method for 
generating the anticontinuous coding sequence for standing- wave phonons from a circle map is described in p7[); in 
this case it is also useful to consider the periodic repetition of sites with ip n = ±a as defects of the zero-amplitude 
state. 

The limit case of one isolated zero-amplitude defect, which is a discrete counterpart of the dark-soliton solution 
of the continuum NLS equation, was investigated in Jic|| . The linear stability analysis of this mode showed that, 
although it is stable close to the anticontinuous limit, it suffers a bifurcation for C/\a\ 2 — C c xs 0.0765 where two 
pairs of eigenvalues of the eigenvalue problem ( pp| ) go out in the complex plane. The resulting oscillatory instability 
occurs due to a resonance between a mode localized around the defect (the defect pinning mode) and linear radiation 
modes. It was shown that for finite systems, the mode recovers its stability above some upper critical value of C/|a| 2 
(since the wavelength of the resonating linear modes becomes larger than the system size) ; however this critical value 
increases with system size so that in the limit of an infinite system, the instability persists for all C/\a\ 2 > C c but 
with a growth-rate that decreases in an exponential- like fashion when approaching the continuum limit C/\a\ 2 — > oo. 
This instability was shown to result in the defect becoming mobile (in NLS terms, the stationary 'black' soliton with 
zero minimum intensity transforms into a moving 'grey' soliton with non-zero minimum intensity) and radiation being 
emitted. In terms of the phase dynamics, this describes a moving, slowly spreading phase kink. 

The instability scenario for the standing-wave phonons is basically the same as for the isolated defect, with the 
essential difference that the localized pinning modes associated with the individual, periodically repeated defects now 
will form a continuous 'defect band'. In general, the Krein signature Q of this defect band is opposite to that of the 
bands associated with the non-zero amplitude sites ]27|] (for < q < ir/2 there are generally several defect bands, 
but they will have the same Krein signature), and as a consequence resonances between the bands will occur if the 
linear coupling C is large enough, giving rise to similar oscillatory instabilities as described above (details are given 
in p7[]). We remark that earlier analysis II]] of standing waves in nonlinear lattices, based on a quasi-continuum 
approximation, did not reveal these instabilities since their origin is the discrete nature of the lattice. 

Let us now return to the main objective of this section, namely to study the effect of the oscillatory standing- 
wave instabilities on the phonobreathers. We find that the families of phonobreathers which are stable close to 
the anticontinuous limit and whose tails approach harmonic standing waves in the small-amplitude limit can be 
constructed from anticontinuous standing- wave solutions at C — 0, placing the breather at a zero- amplitude site of 
the phonon and adjusting it so that the resulting solution is either symmetric or antisymmetric around the breather 
site. Denoting the anticontinuous breather amplitude by b (the phase of the breather site is unimportant when \b\ =/= \a\ 
fj"lf), this yields the following possibilities: 

(i) For q > tt/2 the antisymmetric anticontinuous solution (here q — 2tt/3) 

{^n(0)} = {••• — a, a, 0, —a, a, 0, —a, a, b, —a, a, 0, —a, a, 0...}, (56) 
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with the asymptotic behaviour ip n (0) ~ asm(qn),n — > ±00 in the continuum limit C/|a| 2 — > 00. (Note that this 
solution is antisymmetric only at C = 0.) 

(ii) From (^) we can construct a symmetric solution for g > tt/2 by introducing an additional phase shift of tt at one 
side of the breather site, giving for q = 2tt/3 

{^An(O)} = {••• — a, a, 0, —a, a, 0, —a, a, 6, a, —a, 0, a, —a, 0...}, (57) 

with the asymptotic behaviour ip n (0) ~ acos(qn ± n/2),n — ► ±00 in the continuum limit. 

(iii) For q = Mir/N < tt/2, A" even, the number of consecutive zero-amplitude sites is odd and the antisymmetric 
anticontinuous solution is for e.g. q = ir/4 

{*pn{0)} = {...0, 0, 0, a, 0, 0, 0, -a, 0, 0, 0, a, 0, b, 0, -a, 0, 0, 0, a, 0, 0, 0, -a...} (58) 

behaving as_0„(O) ~ acos(gn), n — > ±00 in the continuum limit. 

(iv) From ( p8[ ) the symmetric solution for q = Mtt/N < n/2, N even, is constructed by a phase shift as above, giving 
for q — 7r/4 

{^„(0)} = {...0, 0, 0, a, 0, 0, 0, -a, 0, 0, 0, a, 0, b, 0, a, 0, 0, 0, -a, 0, 0, 0, a...} (59) 
with the asymptotic behaviour ip n (0) ~ — asin(gn ± tt/2), n — > ±00 in the continuum limit. 

(v) For g = Mtt/N < tt/2, N odd, the number of consecutive zero-amplitude sites is even and we must add an extra 
site to obtain the antisymmetric anticontinuous solution, which for e.g. q = tt/3 becomes 

{V>n(0)} = {...0, 0, a, 0, 0, -a, 0, 0, a, 0, b, 0, -a, 0, 0, a, 0, 0, -a, 0, 0, a...} (60) 

behaving as Vvi(0) ~ acos(q(n + 1/2) ± q/2),n — > ±00 in the continuum limit. (A solution with similar properties is 
obtained by instead removing one zero-amplitude site; however its symmetric counterpart is always unstable.) 

(vi) The symmetric counterpart of (^0|) for q — Mtt/N < tt/2, N odd, is constructed by a phase shift as above, giving 
for q = tt/3 

{^„(0)} = {...0, 0, a, 0, 0, -a, 0, 0, a, 0, b, 0, a, 0, 0, -a, 0, 0, a, 0, 0, -a...} (61) 

with the asymptotic behaviour VVt(0) ~ — asm(q(n + 1/2) ± (g + 7r)/2), n — > ±00 in the continuum limit. 

A typical example on the time-evolution for an initially very weakly perturbed phonobreather with phonon wave 
vector q > tt/2 and phonon amplitude small but non-negligible compared to the breather amplitude is illustrated in 
Fig. H (The example in the figure belongs to type (i), but similar dynamics is observed also for the spatially symmetric 
states of type (ii).) We can clearly distinguish two different steps leading to the final breather destruction. The first 
step is the linear oscillatory instability described above, which leads to the generation of new internal frequencies 
of the breather, and to the movement of the defect sites in a similar way as for the case of an isolated defect. In 
the second step, the moving defects start interacting, and a close inspection of Fig. ^ (a) shows that neighboring 
defects tend to merge and create regions of accumulated phase fluctuations travelling around in the lattice. These 
will interact with the breather, and apparently cause its decay. When the breather has decayed sufficiently to have an 
excitable pinning mode, it will start to move in the lattice but with rapidly decreasing amplitude, and it will finally be 
destroyed. We have at present no complete understanding for the mechanism by which the interaction of the breather 
with the moving 'phase kinks' cause its decay, but we remark that a similar scenario was observed when adding to the 
DNLS equation an external, parametric white noise term p2| . In that case, the white-noise approximation allowed a 
qualitative understanding of the breather decay as a consequence of phase fluctuations by using a collective coordinate 
approach. 

However, to observe this scenario for breather destruction it is necessary (at least for a finite system) that the 
phonon amplitude is not too small compared to the breather amplitude. If we increase the breather amplitude in 
Fig. U (or decrease the phonon amplitude) sufficiently, we find that although the oscillatory instability develops, the 
fluctuations created in the second step will be too weak to cause the breather to decay, and it will live seemingly 
forever as a 'chaotic phonobreather'. The absence of decay for small perturbations can be viewed as a consequence 
of the fact that the single-site DNLS breather is nonlinearly (Lyapunov) stable for norm-conserving perturbations, in 
the sense that \ip n (t)\ remains arbitrarily close to the breather for all times if the initial perturbation is small enough 
[pTstl . Thus, it is clear that for finite systems, the breather cannot be destroyed unless the phonon amplitude exceeds 
some critical value, while nothing can be said about the infinite system since any infinitely extended phonon obviously 
has an infinite norm. 

With < q < tt/2, the first step resulting from the oscillatory instability occurs in a similar way as for q > tt/2: 
the breather acquires new frequencies and the small-amplitude sites of the phonon start moving. By instead making 
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the interpretation that the sites with non-zero amplitude at the anticontinuous limit are defects in the zero-amplitude 
state, their movement can be seen as a consequence of the repulsive interaction between spatially separated, small- 
amplitude breathers with opposite phases observed e.g. in ]l7p . We also observe, similarly as for q > ir/2, the merging 
of neighboring defects, but in this case their interaction with the breather will not lead to breather decay, but rather to 
breather growth if the original amplitude of the breather is not too large. A typical example is illustrated in Fig. ^. It 
can be seen by a careful inspection of Fig. (a) how the merging of small-amplitude sites results in localized humps of 
larger amplitude reminiscent of small-amplitude moving breathers travelling around in the lattice. The interaction of 
these humps with the original breather leads to growth of the latter in a similar way as observed in p9[-pi[. However, 
this growth stops when the breather amplitude has reached a critical value which is close to (but apparently smaller 
than) that corresponding to the limit value A = 4C for small-amplitude perturbations found in Sec. |y| (the latter 
corresponds to |VVt | 2 ~ 5.65). The final state appears also here to be a 'chaotic phonobreather'; we have followed 
the time evolution of this kind of state for times up to 10 6 without seeing any signs of decay. Also, if the initial 
breather frequency is chosen above the critical value A = 4C, we typically do not observe breather growth; instead the 
mean value of the chaotic amplitude oscillations resulting from the oscillatory instability remains close to the initial 
amplitude. 



VI. CONCLUDING REMARKS 



Investigating the interaction between discrete nonlinear Schrodinger breathers and small perturbations, we have 
found firstly that exciting an internal mode of the breather always leads to a slow energy transfer to the breather, 
i.e., to breather growth. Furthermore, we found that a DNLS-breather can pump energy from a small-amplitude 
standing- wave phonon, provided that the phonon wave vector is smaller than the critical value q c given by Eq. Q). 
In both cases, the mechanism for breather growth involves the higher-order generation of radiating modes. Since this 
mechanism disappears at the threshold value A = 4C of the breather frequency, it is impossible for a breather to grow 
beyond this value with the type of small-amplitude perturbations considered here. To analyze the interaction between 
breathers and standing-wave phonons of small but non-negligible amplitude, we considered the long-time evolution 
of weakly perturbed exact phonobreather solutions. The instabilities of these, originating in oscillatory instabilities 
of the nonlinear phonons, where shown to lead to propagating inhomogeneities whose interaction with the breather 
provided a mechanism for breather decay and destruction (when the phonon wave vector q > tt/2) or growth (q < n/2 
andA<4C). 

As was mentioned already in the introduction, the existence of the two conserved quantities (|l|) and (^) makes the 
DNLS equation non-generic among nonlinear lattice equations, and it is therefore necessary to investigate to what 
extent the results obtained in this paper apply also for Klein-Gordon and FPU lattices. We plan to address these 
questions in a forthcoming publication, but let us stress already here that the perturbational approach used here for 
the DNLS equation needs to be modified to account for the fact that generically, the dynamics of the breather involves 
also higher harmonics of its fundamental frequency. Moreover, the approach used in Sees. Ill C and IV based on the 
conservation laws cannot be directly applied in the absence of a second conserved quantity. However, in view of the 
wide applicability of the DNLS equation (and in particular its appearance as a limit case of general lattice equations 
as mentioned in the introduction), we believe that the mechanisms for breather growth and destruction described in 
this paper are essential ingredients also for the corresponding processes in general lattice models. 
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FIG. 1. Variation of internal mode frequencies versus breather frequency for the spatially symmetric (solid line) resp. anti- 
symmetric (dashed-dotted line) internal modes of the single-site breather. Dashed straight line shows the lower band edge of 
the phonon band. 
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FIG. 2. Time-evolution of a breather with an initial perturbation in the direction of the breathing mode. Parameters are 
A = 0.5, uj p ~ 0.47, and C = 1. (a) shows the time-evolution of the central-site intensity IV'nol 2 , (b) shows its time-average 
(iV'nol 2 ) calculated using Eq. (|29[), while (c) shows the instantaneous shift of breather frequency A(t) — Ao. The solid line in 
the main figure in (c) is a fit using Eq. @ with a(0) = 0.082, C\ = 0.9078, C 2 = 0.069, and 7 = 0.067. 
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FIG. 3. Time-evolution of a breather with an initial perturbation in the direction of the pinning mode. Parameters are 
A = 0.45, uip ~ 0.197, and C = 1. (a) shows the time-evolution of |i/>„ _i| 2 , where no is the central site of the breather, (b) 
(main figure) shows the instantaneous shift of breather frequency A(t) — Aq, and inset in (b) shows the time-average (\ip„ | 2 ) . 
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Fig. 4(a), Johansson, Phys. Rev. E 
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FIG. 4. Time-evolution of a breather with an initial perturbation in the direction of an extended eigenmode (spatially 
symmetric) with q < q c . Parameters are A = 1.0, ui p « 2.31 (q ~ 1-22), a ~ 0.0383, and C = 1. (a) shows the time-evolution of 
the central-site intensity \ip„ | 2 , (b) (main figure) shows the total norm contained in a region of 120 sites around the breather, 
and inset in (b) shows the time-average {\^>n \ 2 ) t - 
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Fig. 5(c), Johansson, Phys. Rev. E time Fig. 5(d), Johansson, Phys. Rev. E 



FIG. 5. Time-evolution of a breather with initial perturbations in the direction of extended eigenmodes with q > q c . 
Parameters are A = 1.0, C = 1, and in (a), (b) oj p « 4.93 (q w 2.88) and a « 0.0644, resp. in (c), (d) uj p « 2.72 (g w 1.43) and 
a w 0.0383. (a), (c) show the time-evolution of the central-site intensity I'i/'nol 2 , while (b), (d) show its time-average. 
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FIG. 6. Time-evolution of the phonobreather (^) with q = 2n/3 for a system of 120 sites (periodic boundary conditions), 
perturbed only by the numerical truncation errors. The breather frequency is A& = 1.55, phonon frequency A p h = —2.95 
(phonon amplitude a ~ 0.2), and C = 1. (a) shows | "0"- (*) | 2 t while (b) (main figure) shows the intensity of the breather central 
site l^nol 2 - Inset in (b) shows the inverse participation number R — Af~ 2 ^2 |^n| 4 , which gives a qualitative measure of the 
degree of localization. 
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Fig. 7(a), Johansson, Phys. Rev. E 
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FIG. 7. Time-evolution of the phonobreather (|5£|) with q = 7r/4 for a system of 120 sites (periodic boundary conditions), 
perturbed only by the numerical truncation errors. The breather frequency is A& = 2.2, phonon frequency A p h = —0.5 (phonon 
amplitude a ?s 0.3), and C = 1. (a) shows \ip n (t)\ 2 , while (b) shows the intensity of the breather central site |^n | 2 - 
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